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$-h , Abstract. A new universality class of quantum criticality emerging in itinerant 

& \ electron systems with strong local electron correlations is discussed. The quantum 

criticality of a Cc- or Yb-valence transition gives us a unified explanation for 
unconventional criticality commonly observed in heavy fermion metals such as 
YbRli2Si2 and /3-YbAlB4, YbCu5_ x Al x , and Cclrlns. The key origin is due to the 
locality of the critical valence fluctuation mode emerging near the quantum critical 
end point of the first-order valence transition, which is caused by strong electron 
' correlations for f electrons. Wider relevance of this new criticality and important 

future measurements to uncover its origin are also discussed. 
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1. Introduction 

Quantum critical phenomena in strongly correlated electron systems, which do not 
follow conventional quantum critical phenomena for spin fluctuations [1, 2, 3, 4], have 
attracted much attention. At low temperatures in the paramagnetic-metal phase in 
YbRh 2 Si2 [5, 6] and YbRh 2 (Sio.9sGeo.o5)2 [7, 8], uniform magnetic susceptibility is 
enhanced as x(T) ~ S p ec jg c heat coefficient shows a logarithmic divergence 

C/T ~ — logT for 0.3 K < T < 10 K, giving rise to large Wilson ratio i?w ~ 17.5 
at T = 90 niK [8] much larger than the conventional value of i?w ~ 2 which is 
expected in systems with locally strong correlation. In a wide temperature range for 
20 mK < T < 10 K, T-linear resistivity appears [7]. 

To understand the unconventional criticality, some theoretical efforts have been 
made. Local criticality theory has been intensively discussed [9, 10] and theory of 
tri critical point has been proposed [11]. However, the mechanism and origin of the 
unconventional criticality do not seem to be fully clarified. Furthermore, new aspects 
which seem to be related to the origin of the unconventional criticality have been revealed 
from experimental side. 

Recently, very similar unconventional criticality has been observed in the heavy 
fermion metal /3-YbAlB 4 [12, 13, 14, 15]. Uniform magnetic susceptibility and specific 
heat show x{T) ~ T~ 0,5 and C/T ~ — logT at least for 0.3 K < T < 2 K, respectively, 
giving rise to the enhanced Wilson ratio i?w ~ 25 at T = 0.4 K [14, 15], and T- 
linear resistivity appears in a wide temperature region for 1 K < T < 4 K [12]. Hard 
X-ray photoemission measurement has revealed that the valence of Yb is +2.75 at 
T = 20 K, indicating that /3-YbAlB 4 is an intermediate valence material, which suggests 
the importance of Yb- valence fluctuations [16]. 

The isostructural first-order valence transition occurs in the temperature-pressure 
(P) phase diagram in Ce metal well known as 7-a transition where the critical end point 
(CEP) is located at T CEP ~ 600 K and P ~ 2 GPa [17]. At x ~ 0.1 in Ceo.9_xTho.1La,,, 
Tcep is suppressed to be close to T = K. Although Tcep seems to be not exactly 
zero but finite, at low temperatures T-linear resistivity appears and uniform magnetic 
susceptibility is enhanced. The Wilson ratio i?w ~ 3 is not extremely large, but is larger 
than 2 [18], which is expected to be realized in the system with strong local electron 
correlations. These suggest the importance of critical Ce-valence fluctuations for the 
unconventional criticality. 

In this paper, we discuss that quantum criticality of Yb- and also Ce-valence 
fluctuation is a key origin of these unconventional criticality. We outline a theoretical 
framework for quantum critical phenomena of the valence transition. Then, we show 
that unconventional criticality commonly observed not only in YbRh 2 Si2 and /3-YbAlB 4 
but also in other materials such as YbCus-^ A\ x and Celrlns can be naturally explained 
by quantum valence criticality. We show that fundamental properties of the materials 
are understood coherently from this viewpoint. We discuss that quantum valence 
criticality offers a new class of universality in itinerant electron systems with strong 
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local electron correlations. 

2. Theory of quantum critical phenomena of valence fluctuations 

2.1. Model 

We discuss electronic states of Ce- and Yb-based heavy fermion systems on the basis 
of a generalized periodic Anderson model with inter-orbital Coulomb repulsion as the 
simplest minimal model: 

%pam = e ^ c Lcka + £fE n . + E (VkfL c ka + H.C.) + U n^n 1 ^ 

kcr i<T kcr i 

+ U ic J2nlnl„ (1) 

icTCT 1 

where fi a and Cj CT are the annihilation operator of the f electron and the conduction 
electron at the i-ih site with a spin a, respectively. Here, nf a (a = f , c) are defined 
by n\ a = f- a fi a and n\ a = c\ a Ci a , respectively, is the energy band of conduction 
electrons. S{ is the f level. Vk is the c-f hybridization. U is the onsite Coulomb repulsion 
between f electrons. Uf c is the inter-orbital Coulomb repulsion, which is considered to 
be important in causing the first-order valence transition (FOVT). Since Ce 3+ is a Ai 1 
configuration and Ce 4+ is a 4f° configuration, eq. (1) describes the electronic states for 
Ce-based heavy fermion systems. On the other hand, since Yb 3+ is a 4f 13 configuration 
and Yb 2+ is a 4f 14 configuration, which is a closed shell, if we take the hole picture 
instead of the electron picture, eq. (1) describes electronic states for Yb-based heavy 
fermion systems. 

The model Hamiltonian (1) may be regarded as a generalization of the Falicov- 
Kimball model [19] in which the c-f hybridization is missing, and also a specified 
version of those discussed by Varma at early stage of research on valence fluctuation 
phenomena [20]. A variety of previous works based on Hamiltonian (1) were briefly 
reviewed in refs. [21, 22]. Similar effects of the d-p Coulomb interaction in the so-called 
d-p model have been investigated as a possible charge-transfer fluctuation mechanism 
of the high T c superconductor [23, 24]. 

2.2. Phase diagram 

To clarify the fundamental properties of the model "Hpam, intensive efforts have been 
made: Mean- field theories were applied to %-pam in the d = 3 [25], d = 2 [26], and 
d = 1 [27] systems with d being a spatial dimension. Numerical calculations by the 
DMRG in the d = 1 system [27] and the DMFT in the d = oo system [28] were also 
performed. The phase diagram of "Hpam for a typical set of parameters of Ce- and Yb- 
based heavy fermion systems is shown schematically in Fig. 1. The FOVT surface (dark 
surface) extends to the sharp valence crossover surface (light surface). The critical end 
line is formed between them, which is the edge of the first-order transition surface. The 
quantum critical end point (QCEP) of the FOVT is the point at which the critical end 
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Figure 1. (color online) T-e{-Uf c phase diagram of Hpam- The first-order valence 
transition surface (dark surface) extends to the valence-crossover surface (light surface). 
The critical end line (solid line) touches T = K at which quantum critical end 
point (QCEP) is realized. Superconductivity emerges near QCEP as illustrated by the 
hatched area (see text). Note that the location of the QCEP depends on the strength 
of the c-f hybridization in (1). When the c-f hybridization becomes smaller, the QCEP 
shifts to the smaller {/ fc side. 



line touches the £f-£/f c plane of the ground state. At the critical end line as well as the 
QCEP, valence fluctuation susceptibility 

Xv = "&7 (2) 

diverges, i.e., Xv = oo. Here, rif = Y^ia{ n \ a ) /N s with N s being the number of sites. In the 
small-£f limit, f electron number per site rif becomes nt = l, which is called the Kondo 
regime, corresponding to Ce +3 (Af 1 electron per site) or Yb +3 (4f x hole per site) state. As 
£f increases, rif decreases, and at the FOVT surface, rif shows a discontinuous decrease, 
while at the valence crossover surface, rif shows a continuous but sharp decrease with 
enhanced valence fluctuation. When £f further decreases, rif further decreases, which is 
called mixed the valence regime and corresponds to the Ce +3+5 or Yb 4 " 3 "* 5 state (5 is 
a positive number smaller than 1) i.e., intermediate valence state. Note that rif at the 
FOVT surface and valence- crossover surface depends on the details of model parameters 
in eq. (1) such as momentum dependence of c-f hybridization Vk and band structures 
£k- Namely, at the FOVT as well as valence- crossover surface, rif can be close to 1 and 
also can be 1 — 5. Hence, note that the terms "Kondo" and "mixed valence" in Fig. 1 
merely represent a relatively larger-rif state and relatively smaller-rif state. 

The critical valence fluctuation arising from the QCEP as well as the critical end 
line in Fig. 1 causes several anomalies: Valence fluctuation mediated superconductivity 
appears in the Kondo regime near the QCEP [21, 27] as shown by the hatched area 
at T = K in Fig. 1. Near the QCEP, residual resistivity is enhanced [29] and T- 
linear resistivity appears [30]. Heavy fermion metals CeCu 2 Ge 2 [31], CeCu 2 Si 2 [32], 
and CeCu 2 (Ge x Sii_ x ) 2 [33], where these anomalies were observed, are considered to be 
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located in the valence crossover regime with the superconducting phase in Fig. 1. Since 
inter-orbital repulsion Uf c is considered to have the inter-site character in Ce- and Yb- 
based compounds, most of the compounds seem to have moderate values of U{ c , which 
are located in the valence crossover regime in Fig. 1. For details, readers can refer to 
refs. [22, 34]. 

Since the critical end line extends from the QCEP to the negative temperature 
region in Fig. 1, the CEP can be regarded to be located at negative temperatures for 
Uf c < Ufc CEP with Ufc CEP being Uf c at the QCEP. Namely, the valence crossover line 
for T > (on the light surface in Fig. 1) is extended from the virtual CEP of the 
FOVT at negative temperatures (see Fig. 2(c)). This CEP at negative temperature is 
an analogous concept for the negative Curie or Neel temperature of nearly ferromagnetic 
or antiferromagnetic materials. We note that in the valence crossover regime in Fig. 1, 
the CEP at negative temperatures is moved up toward the T > direction by applying a 
magnetic field, giving rise to an emergence of the QCEP [35] . Namely, it has been shown 
that the CEP is induced to emerge from the negative temperature side to the positive 
temperature side by applying a magnetic field in the valence- crossover regime [34, 36]. 
Correspondence between the phase diagram in Fig. 1 and the T-P or T-H phase diagram 
in Ce- and Yb-based systems is explained in detail in ref . [35] . 

2.3. Quantum valence criticality 

In model (1), the on-site Coulomb repulsion U for f electrons is the strongest interaction. 
Therefore, we should first take into account the effects of the U term. Then, the 
mode-coupling theory for critical valence fluctuations caused by the U& term should 
be constructed. To construct such a theoretical framework, we employ the large-iV 
expansion scheme. Namely, a =t, 4 is generalized to m = 1---N in eq. (1). By 
imposing a constraint on f-electron number by introducing the Lagrange multiplier Aj 
as ~~ J2i \ fem flmfim + Nb\bi — with 6j being a slave-boson operator, we first take 
the U = oo limit in eq. (1). Based on the saddle-point solution for U = oo, the action 
of the system is obtained as a perturbation expansion with respect to the Uf c term 
by introducing the identity applied by a Stratonovich-Hubbard transformation e _s " = 



J%exp[E im / ^r{-|^,(r) 2 + i%{c m f im - / im 4>im(r)}] to S> = jg drC'(r) 
with £ = -t/fc Y.im(ji c im +n\ m ) /2 + U ic E iW n\ m n c im ,/N. Namely, the action S = S + S' 



where Sq is the action for the terms without Uf c in (1) is given by 



S M = oE ^2{q)Vm{q)¥m{-q) 



in q 
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where fi 2 (g) « U {c [l - 2£/" fc Xo cc (g)/iV] , fi 3 (0i,02,03) = 0(C/|), and Ot(0i, 02, 03, 04) = 
0(Uf c ) with g = (q, iwi) and = 2/7rT. Here, Xo cc (0) is given by Xo° c (0) = 
-TE k , n ^(k + q,i£ n + z^)G^ c (k,ze„)/iV s with e„ = (2n + l)vrT. Here, G*(k,ie n ) 
and GQ C (k, ie n ) are the Green functions of f and conduction electrons for the saddle 
point solution for U = 00, respectively [37], and N s is the number of lattice sites. 

Since long wavelength |q| <C q c around q = and low frequency \u\ <C u c 
regions play dominant role in critical phenomena with q c and u c being cutoffs for for 
momentum and frequency, respectively, fi, are expanded for q and u around (0, 0) as 



^2(0) ~ V + M 2 + C q \uJi\, n 3 (0i, 02, 03) ~ v 3y jT/N s , and Q 4 (0i, 02, 03, 04) ~ v±T/N s . Let 
us here apply the Hertz's renormalization-group procedure [3] to S[<p\: (a) Integrating 
out high momentum and frequency parts for q c /s < q < q c and u c /s z < uj < oj c , 
respectively, with s being a dimensionless scaling parameter (s > 1) and z the dynamical 
exponent, (b) Scaling of q and u by q' = sq and uj' = s z u. (c) Re-scaling of <p by 
ip'(q',uj') = s a ip(q' /s,u' /s). Then, we dertermined that to make the Gaussian term in 
Eq. (3) scale invariant, a must satisfy a = —(d+z+2)/2 with d spatial dimension and the 
dynamical exponent z = 3. The renormalization-group equations for coupling constants 
Vj are derived as ^ = [6 - (d + z)} v 3 + 0(v%), and ^ = [4 - (d + z)] v A + 0(vf), for 
cubic and quadratic terms, respectively. By solving these equations, it is shown that 
higher order terms than the Gaussian term are irrelevant 

lim Vj(s) = for 7 > 3 (4) 

for d + z > 6. For the case of d = 3 and z — 3, it is shown that the cubic term is 
marginally irrelevant [22]. It is noted that in d = 3 systems for realistic heavy fermion 
systems, eq (4) is considered to hold as discussed in the end of section 2.3. 

We have found that almost dispersionless critical valence fluctuation modes appear 
near q = not only for deep £f, i.e., in the Kondo regime, but also for shallow Ef, i.e., 
in the mixed valence regime, because of strong on-site Coulomb repulsion for f electrons 
in eq. (1) [37]. This "almost local" nature causes the extremely small coefficient A in 
dynamical valence susceptibility defined by the inverse of the coefficient of the Gaussian 
term f2 2 (0) in the action 

Xv(q,m) = : , 2- n 1 r , (5) 
n + Aq 2 + C q \ui\ 

where C q = C/max{g, Z;" 1 } with Zj being the mean free path of impurity scattering [38] 
and bj\ is the Boson Matsubara frequency ui = 2nlT. This yields C q = C/q in the clean 
system and C q = C/k in the dirty system affected by impurity scattering, giving rise to 
the dynamical exponent z = 3 and z = 2, respectively. 

The emergence of the weak-g dependence in the critical valence fluctuation is 
analyzed, as follows [22]: Since Xvio) — ^2~ 1 (0), ^ ne ^-dependence in eq. (5) appears 
through Gq (k + q) in Xo° c (0)- Near q — 0, Xo c °(0, 0) is expanded as 

Xo CC (0,O) = Xo CC (O,O) + sf^— ) q\ (6) 

VIA* — £ fl/ 
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where S includes the effect of the f-electron self-energy for U in eq. (1). Since the 
f-electron self-energy has almost no q dependence in heavy electron systems, the q 
dependence of the f-electron propagator Gq (k + q) comes from the hybridization V 
with conduction electrons with the dispersion £k+ q , as seen in the coefficient of the 
q 2 term in eq. (6). Hence, the reduction of the coefficient A in eq. (5) is caused by 
two factors. One is due to the smallness of {V/\fi — E{\) 2 . In typical heavy electron 
systems, this factor is smaller than 10" 1 . The other is the reduction of the coefficient 
S, which is suppressed by the effects of the on-site electron correlations U in eq. (1). 
Numerical evaluations of Xo cc (<?'0) based on the saddle point solution for U = oo in 
eq. (1) show that extremely small S appears not only in the Kondo regime, but also in 
the mixed- valence regime [37], indicating that the reduction by S plays a major role. 
These multiple reductions are the reason why extremely small coefficient A appears in 



The extremely small A in eq. (5) makes the characteristic temperature for critical 
valence fluctuations 

extremely small. Here, q^ is a momentum at the Brillouin zone boundary. Hence, even 
at low enough temperature, lower than the effective Fermi temperature of the system, 
i.e., the so-called Kondo temperature, T <C Tk, the temperature scaled by To can be 
very large: t = T/T ^> 1. This is the main reason why unconventional criticality 
emerges at "low" temperatures, which will be explained below. 

By optimizing the action (3) derived from Upam, taking account of the 
mode coupling effect for critical valence fluctuations, we obtain a self-consistent 
renormalization (SCR) equation for critical valence fluctuations: 



where y = r]/(Aq^), x = q/qs, x c = q c /qB with q c being a momentum cutoff, and y 
parameterizes a distance from the criticality and y\ is a dimensionless mode-coupling 
constant of 0(1). The solution of Eq. (8) is quite different from that of ordinary SCR 
equations for spin fluctutions [1] because of the extremely small A in eq. (5). 

In the y 3> t limit, the analytic expression of Xv(0,0) can be obtained, which 
is Xv(0,0) = y^ 1 ~ t~ 2 / 3 for both the clean (z = 3) system and dirty (z = 2) 
system. As shown in ref. [36], uniform magnetic susceptibility diverges at the QCEP 
where valence fluctuation diverges. This is due to the fact that Xv(q, iwi) and the 
dynamical f-spin susceptibility %(q, iui) = Jq dr(T T S^(q, r)Sf(— q, r))e lu)lT have the 
common structure near the QCEP [35]. The spin-lattice relaxation rate is given 
by (TiT)" 1 = 2(7 2 /5( f/ u B ) 2 2 q |D q | 2 Imx R (q,w )/w , where x R (q,^o) is the retarded 
dynamical f-spin susceptibility with the nuclear resonance frequency ojq, 7 n is the 
gyromagnetic ratio of nuclear spin, is a Lande's g factor for f electrons, and D q 
is the hyperfine-coupling constant. Then, we have (TiT) -1 w y^ 1 near the QCEP. 



eq. (5). 




(8) 
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Table 1. Resistivity, specific-heat coefficient, uniform magnetic susceptibility, and 
spin-lattice relaxation rate in the vicinity of the QCEP of first-order valence transition. 
Exponent a takes the value for 0.5 < a < 0.7 depending on the flatness of the critical 
valence fluctuation mode and measured T range (see text). 

p C/T X (TiT)- 1 

T -logT T~ a T~ a 



Namely, the uniform magnetic susceptibility and spin-lattice relaxation rate are shown 
to have the same temperature dependence: x(t) ~ t~ 2 ^ 3 and (TiT)" 1 ~ t~ 2 / 3 for y ^> t. 
When T is decreased down to T ~ T , y in eq. (8) is evaluated as y ~ t ' 5 by the 
least square fit of the numerical solution of eq. (8). Hence, depending on the flatness 
of critical valence fluctuation mode and measured temperature range, x{T) ~ t~ a and 
(TiT)^ 1 ~ t~ a with 0.5 < a < 0.7 are observed. 

Calculating resistivity p{T) as p{T) oc ^ dujujn{uj)[n{ui) + 1] Jq c dqq 3 lmx^(q, u>) 
with n{uj) = l/(eP u — 1) being the Bose distribution function, and Xv(9i u ) = 
(n + Aq 2 — iCqu)' 1 , a retarded valence susceptibility, we obtain p(t) oc t for the y > t 
region. Here, emergence of T-linear resistivity comes from the high-temperature limit 
of the Bose distribution function, indicating that the system is described as if it is in the 
classical regime, because the system is in the high-T regime in the scaled temperature 
t = T/Tq ^> 1, in spite of T Tk- Hence, this mechanism is essentially the same as 
the emergence of T-linear resistivity in the electron-phonon system where p(T) ~ T 
appears for T > Od/5 with 6d being Debye temperature [39]. This result that the 
almost flat dispersion of the critical valence fluctuation mode causes T-linear resistivity 
implies that dynamical exponent z can be regarded as if z = oo. 

An evaluation of the valence-fluctuation exchange process in the self-energy by 
using eq. (3) gives logarithmic energy dependence in the real part of the self-energy [42] . 
A numerical solution of the self-energy gives logarithmic temperature dependence in 
the specific-heat coefficient C/t ~ — logt even for the temperature regime for T > T 
where the uniform magnetic susceptibility and the resistivity behave as x{t) ~ t~ a with 
0.5 < a < 0.7 and p(t) ~ t, respectively [36]. Anomalous temperature dependence 
of the specific-heat coefficient and the resistivity of our theory shares an aspect 
similar to that based on the marginal Fermi liquid (MFL) assumption for high-T c 
cuprates [40, 41], in which the self-energy is given by oc e (Ins — i|e|). Since uniform 
magnetic susceptibility x(T) shows a power-law divergence, the Wilson ratio Rw, i.e., 
x/(C/T), is significantly enhanced: i?w ^ 2. When the experimentally accessible 
lowest temperature is larger than To, the unconventional criticality dominates all the 
physical quantities down to the lowest temperature. Namely, near the QCEP in Fig 1, 
unconventional criticality summarized in Table 1 appears. 

When the experimentally accessible lowest temperature is smaller than T , y oc t 4 / 3 
(y oc t 3 / 2 ) is realized in the clean (dirty) system for the t = T/T < 1 region. Hence, 
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Figure 2. (color online) Focus on the vicinity of T = K in T-P phase diagram with 
first order valence transition (solid line) terminating at critical end point (CEP) (filled 
circle) from which sharp valence crossover extends (dashed line), (a) Tcep > case, 
(b) T CEP = case, i.e, quantum critical end point (QCEP) is realized, (c) T CEP < 
case. Note that in Ce (Yb)-based systems, the Kondo regime and mixed-valence regime 
correspond to the smaller (larger) P regime and larger (smaller) P regime, respectively. 
This T-P phase diagram corresponds to the cutout of the T-gf plane in Fig. 1 at 
moderate- t/fc region, where the CEP is located at negative temperature (see ref. [35] 
for detail). 



resistivity behaves as p ~ t 5 ^ 3 in the clean system and p ~ t 3 ^ 2 in the dirty system in 
the t -»■ limit (T < T < T K ). 

We note that the present formulation reproduces the previous results for the 
local limit of critical valence fluctuation mode A — > in eq. (5), which was derived 
starting from phenomenologically introduced dynamical valence susceptibility assuming 
z = oo [30]. 

By the third law of thermodynamics, entropy of the system should be zero at the 
ground state. This requires that the FOVT line appears in the direction parallel to the 
temperature axis of the phase diagram. For example, in the case of the T-pressure (P) 
phase diagram, the slope of the FOVT temperature T v is given by 

dT v _ Vk — Vmv ^ 
dP Sk — Smv 

by the Claudius- Clapeyron relation. Here, V and S denote the volume and entropy, 
respectively. The subscript K and MV are abbreviations for the Kondo and mixed- 
valence regimes, respectively, which represent the relatively larger-rif state and smaller- 
rif state, respectively, as remarked in §2.2. By the third law of thermodynamics, the 
entropy should be zero at T = K. Hence, the denominator of eq. (9) should be zero 
at the ground state, giving rise to dT w / dP\r=o = oo. This indicates that the FOVT line 
should appear perpendicularly to the P cLXIS, clS shown in Fig. 2(a). Here, the valence 
crossover line extends from the critical end point of the FOVT line. Hence, the valence- 
crossover line extended from the QCEP appears perpendicularly to the P axis, as shown 
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in Fig. 2(b). Then, when T is decreased to approach the QCEP, we usually follow the 
valence crossover line shown by the dashed line in Fig. 2(b) at least in the vicinity of 
the QCEP. 

At the valence crossover line, the cubic term in the action (3) derived from "Hpam 
eq. (1) [37] in the vicinity of the QCEP vanishes, which makes the upper critical 
dimension d u = 4 of the system, but not d n = 6. Then, the clean (z = 3) system 
and the dirty (z = 2) system in three spatial dimension (d = 3) are both above the 
upper critical dimension, i.e., d + z > d n . Then, the higher order terms other than the 
Gaussian term are irrelevant in the action, which makes the fixed point Gaussian. This 
guarantees the validity of eq. (8) derived by constructing the best Gaussian in the action 
for "Hpam- Note that this conclusion holds even when P is replaced by the other control 
parameters such as a magnetic field H in Fig. 2. 

2.4- Valence fluctuation mediated superconductivity 

Although almost flat dispersion appears near q = for critical valence fluctuation mode 
as presented above, q dependence of the critical mode, i.e., a gradual decrease, appears 
around q = 2kp with kp being the Fermi wavenumber. This moderate q dependence of 
critical valence fluctuation is shown to induce attractive interaction at nearest-neighbor 
distance in the real space by the Fourie transformation [22]. Since strong Coulomb 
repulsion U in eq. (1) forbids on-site pairing, d-wave pairing is induced by critical valence 
fluctuations for the spin-singlet sector. Indeed, d-wave pairing is shown to appear in the 
Kondo regime near the QCEP in Hpam as shown in Fig. 1 by the slave-boson mean-field 
theory taking account of Gaussian fluctuations in the d = 3 system [21] and the DMRG 
calculation in the d = 1 system [27]. We note that the direction of node in d-wave 
pairing, i.e., d x 2_ y 2 or d xy symmetry, depends on the details of each system such as 
the shape of Fermi surface and the lattice structure. In principle, spin-triplet pairing, 
such as p-wave pairing, is not excluded because it is also included in a manifold of the 
intersite pairings. 

3. Comparison with experiments 

3.1. /3-YbAlB A 

/3-YbAlB4 is a heavy fermion metal with intermediate valence, which is Yb +2 ' 75 at 
T = 20 K [16]. The characteristic energy scale corresponding to the so-called Kondo 
temperature is estimated to be Tk ~ 200 K, which has a large c-f hybridization [12]. 
The uniform magnetic susceptibility and specific heat coefficient show a non- Fermi liquid 
behavior \ ~ T -0 ' 5 and C/T ~ — logT for 0.3 < T < 2 K, respectively, giving rise to 
large Wilson ratio i?w ~ 25 at T = 0.4 K [12, 15]. A T-linear resistivity emerges in 
a wide temperature range for 1 K < T < 4 K and p ~ T 15 for T < 1 K [12]. These 
unconventional criticality summarized in Table 2 are well explained by the quantum 
valence criticality discussed in §2.3 (see Table 1). We note that the nuclear spin- 
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P 


C/T x 


(TiT)- 1 


Rw 


reference 


YbRh 2 Si 2 


T 


-logT T" 06 


rp-0.5 


17.5 


[5, 6, 7, 8, 43] 


/3-YbAlB 4 


rpl.5 y rp * 


-logT T" ' 5 


** 


25 


[12, 13, 14, 15] 



Table 2. Physical quantities of resistivity p (* p ~ T 15 for T < 1 K, p ~ T for 
T > 1 K), specific heat coefficient C/T, uniform magnetic susceptibility x(T), nuclear 
spin-lattice relaxation rate (TiT) -1 i** experiment is desired), and Wilson ratio i?w 
at low temperatures in YbRh.2Si2 and /3-YbAlB4. 

lattice relaxation rate (TiT) -1 ~ T~ a with 0.5 < a < 0.7 is predicted by the theory 
of quantum valence criticality, and the measurement of (TiT)" 1 is highly desirable. 

3.2. YbRh 2 Si 2 

In YbRh 2 Si 2 [5, 6] and YbRh 2 (Si .95Ge .o5)2 [7, 8], unconventional criticality is observed 
as uniform magnetic susceptibility x(T) ~ T~ 6 , specific-heat coefficient C/T ~ — logT 
for 0.3 < T < 10 K [7], giving rise to large Wilson ratio ,R W « 17.5 at T = 90 mK [8], 
T-linear resistivity for 20 mK < T < 10 K [7], and spin-lattice relaxation rate 
(TxT)" 1 ~ f or 50 mK < T < 1 K [43]. These are summarized in Table 2 and are 

well explained by the quantum valence criticality (see Table 1). 

At T = 300 K, the valence of Yb is observed as v = +2.90 [44]. This 
strongly suggests that when T decreases, the Yb valence further decreases to be of 
the order of v ~ +2.8 near T = K. This implies that YbRh 2 Si 2 is an intermediate 
valence material [45], similar to /3-YbAlB4 [16]. Recent angle- resolved photoemission 
measurement at T = 15 K has revealed that the c-f hybridized band exists near the 
Fermi level in YbRh 2 Si 2 [46]. This clearly indicates that the c-f hybridization exists at 
low temperatures in YbRh 2 Si 2 for H = 0. 

In the T-H phase diagram of YbRh 2 Si 2 , a characteristic temperature T*(H) 
appears, as illustrated by a dashed line in Fig. 3(a), where physical quantities show 
crossover behavior [6]. The origin of T*(H) is considered to be closely related to 
the origin of unconventional criticality in YbRh 2 Si 2 . Here, we note that a field- 
induced QCEP of the FOVT can explain the emergence of the T*(H) line. Recently, 
magnetic field dependence of the QCEP has been theoretically studied on the basis of 
"Hpam (1) [25, 36]. It has revealed that the QCEP extends to the smaller-£/f c region in 
Fig. 1 as magnetic field increases, which makes a sharp valence crossover line T*(H) 
appear in the T-H phase diagram, as shown in Fig. 3(a). 

Recently, interplay of the QCEP of the FOVT and magnetism has been theoretically 
studied in "Hpam [26, 35]. For a realistic strength of c-f hybridization |Vk| in eq. (1), 
a coincidence of magnetic transition temperature T mag and sharp valence crossover 
temperature T* occurs, as shown in Fig. 3(a) [26]. When the hybridization decreases, 
H c tends to be larger than H* in Fig. 3(a). Then, enhanced valence fluctuations at 
H = H* suppress the magnetic order, which still makes the coincidence of H c w H*. 
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Figure 3. (color online) (a) Schematic T-H phase diagram in the case of small to 
moderate strength of c-f hybridization in %pam (see text). T mag is the magnetic 
transition temperature and T* is valence crossover temperature. Energy bands for 
(b) paramagnetic metal phase and (c) AF-ordercd phase in Hpam- In (c), the shaded 
dashed line represents energy band for conduction electrons, fcp represents the small 
Fermi surface for only conduction electrons and the solid dashed line is a guide for the 
eyes (see text). 

Namely, for realistic values of the c-f hybridization |Vk| ranging from the rather small 
to moderate strength, the coincidence of T mag and T* occurs at T = K. This result 
well explains the T-H phase diagram of YbRh 2 Si2- 

We also note that the T-H phase diagram shown in Fig. 3(a) was actually realized 
in YbAuCu 4 , where Cu-NQR frequency change was detected at T*(H) clearly indicating 
Yb- valence crossover [54] . The phase diagram is very similar to the T-H phase diagram 
in YbRh 2 Si 2 [34]. 

When we set a larger |Vk| in "Hpam, the magnetically ordered phase shrinks to be 
H c < H* in Fig. 3(a). Namely, the magnetic quantum critical point (QCP) and the 
valence crossover point are separated. On the other hand, when we set very small c-f 
hybridization, H c exceeds H* to be H* < H c . Namely, crossing of T mag and T* occurs. 

In the paramagnetic phase in YbRh 2 Si2, it was demonstrated that a tiny change in 
Yb valence, about 0.03, can reproduce the measured drastic change in the Hall coefficient 
by the band-structure calculation [47]. Namely, this suggests that the measured Hall- 
coefficient change at T*(H) in the T-H phase diagram of YbRh 2 Si 2 can be explained 
by the Yb- valence crossover. In ref. [47], it was pointed out that a slight shift of the 
f level causes large change in the conductivities because of large density of states for 
4f electrons. Here, it is noted that band-structure calculations for YbRh 2 Si 2 have been 
performed by several groups [48, 49, 50] and a renormalized-band approach has been also 
discussed [51]. Although these methods do not seem to show complete agreements about 
detailed band structures and Fermi-surface shapes, they at least show large density of 
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states for 4f electrons. Hence, it seems important to examine the possibility of whether a 
slight shift of the large density of states causes the drastic change of the Hall coefficient 
under a magnetic field within the band-structure calculations. 

Although the Kondo breakdown scenario has been intensively discussed to explain 
the field dependence of the Hall coefficient [52], it would not be clear enough how 
it is reconciled with the experimental fact that large effective mass lim.T^oC/T ~ 
1.7 J/molK 2 is observed at H = [53] (see also §4). Here, we focus on the quantum phase 
transition between the AF- and paramagnetic- metal phases (see H ~ H c in Fig. 3(a)) 
to address this issue from the viewpoint of the Yb-valence crossover. 

In the paramagnetic metal phase, a large Fermi surface which includes an f-electron 
number is realized, as shown in Fig. 3(b). In the antiferromagnetic (AF) order phase, 
the hybridized bands are folded, as shown in Fig. 3(c). An important point here is that 
because the lowest folded band with one electron number per site shrinks down from 
the Fermi level, the folded band at the Fermi level has an electron number which is the 
same as that for conduction electrons. Namely, the Fermi surface in the folded band is 
the same as the small Fermi surface only made from the conduction electrons, as shown 
by fcp in Fig. 3(c). Hence, a small to large Fermi surface change occurs at the AF to 
paramagnetic transition, even though c-f hybridization is always finite, i.e., (f^c^) ^ 0. 
Note that this is just a matter of number counting, which does not depend on the order 
of magnetic transition, i.e., the first-order or second order transition does not matter. 
Because heavy quasi particles contribute to the formation of the AF-ordered state, as 
shown in Fig. 3(c), mass enhancement occurs in the AF phase. This naturally explains 
the fact that heavy electron mass is observed as C/T ~ 1.7 J/molK 2 in the magnetically 
ordered phase in YbRh 2 Si2 [53]. These aspects cannot be understood in the so-called 
Kondo breakdown scenario [9, 10]. 

To clarify the origin of unconventional criticality in YbRh 2 Si2, we propose to 
perform the Co-NQR measurement in the Co-doped sample: Yb(Co x Rhi_. r )2Si2 [55]. 
If the Co-NQR frequency z/q changes at the characteristic temperature T*(H) in the 
T-H phase diagram, it will be evidence that Yb-valence crossover occurs at T*(H) as 
observed in YbAuCu4 [54]. Measurement of this is highly desirable. 

3.3. YbCu 5 _ x Al x 

In the heavy fermion metal YbCu5_ x Al x , at T = 300 K, a gradual decrease in the Yb 
valence occurs as x decreases from 2.0 to 0.0. However, at T = 10 K, a sharp decrease 
in the Yb valence occurs near x = 1.5 [56], which has been also supported by recent 
detailed measurements of the Yb valence [57, 58]. These results indicate that a sharp 
Yb-valence crossover occurs near x = 1.5 with strong Yb-valence fluctuations. The 
schematic T-x phase diagram of YbCus-^AL is shown in Fig 4. At x — 1.5, as T 
decreases, logarithmic divergence in the specific- heat coefficient C/T ~ — logT and the 
power-law divergence in the uniform magnetic susceptibility x ~ T _2//3 were observed 
in 1997 by E. Bauer et al [56], which well agree with the quantum valence criticality 
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Figure 4. (color online) Schematic phase diagram of YbCus-^Ala; [56, 57, 58]. Near 
x = 1.5 sharp Yb- valence crossover occurs between the Kondo regime and mixed 
valence regime at low temperatures T < 10 K. 



(see Table 1). 

The shape of T*(g), i.e., the control parameter dependence of the valence crossover 
temperature with g being the control parameter such as a magnetic field, pressure and 
chemical substitution, depends on where the material is located in the phase diagram 
in Fig. 1. Detailed comparison of T*(g) between the theory and experiments will be 
interesting future studies. 



34. CeMn 5 

Celrln 5 is a heavy fermion metal which exhibits superconductivity at T c « 0.4 K [59]. 
Interestingly, by applying pressure, spin fluctuations are quickly suppressed, while 
superconducting transition temperature T C (P) increases and has a maximum around 
P ~ 2.5 GPa, which suggests that the origin of superconductivity may be different 
from ordinary spin-fluctuation mechanism [60]. Recently, M. Yashima and Y. Kitaoka 
have observed that In-NQR frequency uq starts to change sharply around P ~ 2.2 GPa 
(namely, z/q — P relation changes its slope for dv^jdP to take a minimum as a function 
of P) [61], as observed in CeCu2Si2 at P ~ 4 GPa [62]. This indicates that Ce- valence 
crossover sharply occurs around the pressure of T C (P) maximum, strongly indicating 
that superconductivity is caused by the Ce-valence fluctuation. 

In the T-H phase diagram, a first-order-transition or sharp-crossover line T*(H) 
appears [63, 64, 65, 66, 67], which is most likely to be FOVT or valence crossover [25]. 
Interestingly, quantum critical phenomena appear near H = H* where the T*(H) line 
touches T = K. As H increases to approach H*, convex curve in the T 1,5 plot of p(T) 
emerges, indicating that T-linear resistivity appears near H = H* = 28 T [66] . Residual 
resistivity po is enhanced around H = H* [67]. Logarithmic divergence in specific heat 
C/T ~ — logT appears near H = H* [64]. We pointed out that field dependence of 
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the QCEP in "Hpam well explains the emergence of the T*(H) line in the T-H phase 
diagram of Celrln 5 [25, 36]. 

We here note that quantum valence criticality shown in Table 1 explains the 
measured criticality near H = H*. Since the QCEP of the FOVT is induced by applying 
the magnetic field [25, 36], quantum valence criticality emerges in high-field region, 
which is consistent with the observations near H ~ H* = 28 T in Celrln 5 [25]. As 
shown in refs. [21, 27], the superconducting phase appears robustly even in the region 
rather far from the QCEP in the phase diagram of "Hpam (see Fig. 1). This is also 
consistent with the recent discovery by Yashima et al, that at H = 0, Celrln 5 seems to 
not be very close to the QCEP, but located in a gradual Ce-valence crossover region. 
Hence, quantum valence criticality based on "Hpam gives a unified explanation for the 
measured anomalies in Celrlns. 

These observations indicate that a new viewpoint of the closeness to the QCEP 
of the FOVT is indispensable for comprehensive understanding of Cell5 systems, in 
addition to the conventional picture of the so-called "Doniach's phase diagram" . Indeed, 
it has been shown theoretically that interplay of the QCEP of the FOVT and magnetism 
gives a unified explanation for several anomalies such as drastic change of Fermi surface 
and non- Fermi liquid transport observed in CeRhIn 5 under pressure [25, 26]. Namely, 
the results discussed for Fig. 3(a) in §3.2 basically hold even when the magnetic field 
is replaced by pressure as a control parameter [26]. For details, readers can refer to 
refs. [26, 35] 

3.5. Discussion 

In CeCu 2 Ge 2 , X-ray absorption measurement under pressure reported that the lattice 
constant shows a sudden jump at P = Py ~ 15 GPa, where the superconducting 
transition temperature has a maximum [68]. However, reexamination of the X-ray 
measurement in CeCu2Ge2 has observed a smooth decrease in the lattice constant under 
pressure but no jump at P = P v [69]. A smooth decrease in the lattice constant around 
P ~ P v ~ 4.5 GPa has been also observed in CeCu2Si2 [69]. We note that a sharp 
Cu-NQR- frequency change has been observed at P ~ P v even without a jump in the 
lattice constant in CeCu2Si2 [62]. This indicates that NQR is more sensitive than the 
lattice-constant measurement, which can detect a tiny change of the electronic state due 
to the Ce- and Yb-valence change with the highest sensitivity. 

Here, let us remind the fact that even in Ce metal [17, 70] and YbInCu4 [71, 72], 
which show isostructural FOVT, the valence change is only about 0.1. In Ce metal, 
isostructural first-order valence transition is known to occur, well known as 7-a 
transition [17]. At T = 300 K, the Ce valence shows a discontinuous jump from Ce +3 ' 03 
to Ce +3,14 at the 7-a transition, as pressure increases. The lattice constant also shows a 
discontinuous jump at the 7-a transition, about 16% shrinkage of the unit-cell volume. 
In YbInCu4, isostructural FOVT is also known to occur [71]. As temperature decreases, 
the Yb valence shows a discontinuous change from Yb +2,97 to Yb +2,84 with a volume 
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expansion, i.e., 0.5% increase in the unit-cell volume, at T v = 42 K [73, 74, 75, 76]. 
Hence, in the valence crossover regime where most of Ce- and Yb-based heavy fermion 
systems are considered to be located, the valence change is expected to be in the order 
of 0.01, even when sharp valence crossover occurs (see Fig. 1). 

Note that the volume change in YbInCu4 is much smaller than that in Ce metal, 
although Ce metal and YbInCu4 both show a jump in Ce and Yb valence of about 
~ 0.1 at the FOVT, respectively. There is a tendency that for volume change in Ce- 
and Yb-based compounds to be much smaller than Ce (or Yb) mono metal. Hence, it 
should be kept in mind that even in the cases where the lattice constant merely shows 
a monotonic change as control parameters, such as pressure and a magnetic field, are 
changed, there is a possibility that NQR frequency change is detectable. 

Recently, the pressure dependence of Ce valence in CeCu2Si2 has been observed by 
X-ray Lm edge absorption spectra [77] . Although the data have error bars of the order 
~ 0.01 and discretely measured pressure data points do not seem to be enough to see 
whether sharp valence crossover in the order of ~ 0.01 occurs around P v ~ 4.5 GPa, a 
possibility of monotonic Ce-valence crossover is not excluded. Since the measurement 
has been performed at T = 14 K, this result may suggest that the QCEP of the FOVT 
is located at much lower temperature, as shown in Fig. 2(c). 

In our studies, we have neglected effects of excited state(s) of crystalline-electric- 
filed (CEF) levels. For the moment, it is not clear whether those CEF states give an 
essential effect for realistic CEF states while an effect of orbital transfer between CEF 
levels (with a model with hybridization with conduction electrons) was discussed by 
Hattori [78]. In any case, such effects certainly deserve more investigations. 

4. Discussion about other theoretical proposals 

In this paper, we have discussed that fundamental properties of the outstanding 
materials such as YbRli2Si2 are explained coherently from the viewpoint of quantum 
valence criticality. So far, to understand the unconventional criticality, several scenarios 
such as the local criticality theory [9, 10, 79] and theory of tricritical point (TCP) [11] 
were proposed. Here, we discuss them within the scope of the present paper. 

The local criticality theory did not fully account for the unconventional criticality 
in YbRli2Si2 shown in Table 2. In ref. [9], it was shown that the magnetic susceptibility 
behaves as x{T) = 1/(6 + AT a ). However, the measured value of a ~ 0.6 was not 
shown explicitly by the theory [9, 79]. The theory concluded that \jT\ = constant, 
which does not reproduce the measurement (TiT)" 1 ~ f or j QW temperatures 

(50 niK < T < 1 K) [43]. Furthermore, the T-linear resistivity and logarithmic-T 
dependence of C/T were not shown by the theory itself [9, 79]. In ref. [10], a scenario of 
Kondo breakdown was just proposed, but the temperature dependences in Table 2 were 
not shown theoretically. In ref. [10], the transition between the conventional SDW-type 
AF metal with reconstructed Fermi surfaces around "hot spots" and the paramagnetic 
metal (PM) with the large Fermi surface is argued. We note here that the "small"- 
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to-large Fermi surface change can occur at the SDW-type AF-PM transition in the 
ground state of the periodic Anderson model for realistic parameters for heavy electron 
systems [26, 80] (see Fig. 3(c) in §3.2). Namely, Fermi surface changes from the "small" 
Fermi surface, which is the same as that of the non-hybridized conduction band, realized 
in the SDW-type AF phase to the large Fermi surface in the PM phase. Furthermore, as 
noted in §3.2, it is not clear how measured large C /T ~ 1.7 J/molK 2 [53] is reconciled 
with switching off the c-f hybridization in the magnetically-ordered phase. Although 
the concept of "dynamical Kondo screening" was proposed to explain the measured 
large mass [79], it seems necessary to show explicitly how large mass can appear in their 
model, i.e., the Kondo lattice model with the RKKY interaction model, for reliable 
argument. It is also noted that in the most of Ce- and Yb-based materials, intermediate 
valences of Ce and Yb are realized near T = K, implying the existence of valence 
fluctuations at low temperatures. 

The TCP theory shows at low temperatures the uniform magnetic susceptibility 
x{T) ~ T~ - 75 , the nuclear spin-lattice relaxation rate (TiT) -1 ~ —-Do logT + DqT~ 75 , 
the resistivity p(T) ~ T L5 , and the Sommerfeld coefficient C/T ~ constat — T 1//2 
(T < 1 K) — log T (T > IK) [11]. Although it may be possible to interpret that 
some agreements with the measurements can be seen for some temperature regions, 
these results are not exactly the same as the measurements shown in Table 2. In 
particular, it should be noted that the TCP theory did not show the T-linear resistivity 
at low temperatures which is the prominent feature in YbRh 2 Si2- Furthermore, it is not 
clear why the characteristic temperature T*(H) appears in the T-H phase diagram of 
YbRh 2 Si2, which is closely related to the origin of the unconventional criticality. 

It is remarked that by the interplay of the magnetically ordered temperature and 
the sharp-valence crossover temperature as a function of control parameters such as 
pressure, magnetic field, or chemical substitution, the first-order magnetic transition 
can be caused, giving rise to the TCP, as shown in ref. [26, 35]. Search for such a 
material is left for future experimental studies. 

5. Summary and Perspective 

By constructing a mode-coupling theory for critical valence fluctuations taking account 
of strong on-site Coulomb repulsion in the periodic Anderson model with inter-orbital 
Coulomb repulsion, we have shown that a new class of universality emerges near the 
QCEP of the FOVT. The quantum valence criticality gives a unified explanation 
for unconventional criticality observed in heavy fermion metals such as /3-YbAlB 4 , 
YbRh 2 Si2, YbCus-zALj;, and Celrln 5 . 

The key origin of emergence of unconventional criticality is clarified to be due to 
the extremely dispersionless critical valence fluctuation mode caused by the strong on- 
site Coulomb repulsion U for f electrons. Observation of the locality of the critical 
valence-fluctuation mode is an important future issue for revealing direct evidence 
of the origin of the unconventional criticality. Resonant X-ray scattering offers a 
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direct observation of q dependence of the valence susceptibility x^il^)- As discussed 
in §2, an almost dispersionless critical valence fluctuation mode is expected to be 
observed near the QCEP of the FOVT. To detect such a tiny width of dispersion, 
high resolution of resonant X-ray scattering measurement is necessary. Since dynamical 
f-spin susceptibility x R {Qi U} ) is considered to have a common structure with the valence 
susceptibility x^(g,u;) [37], neutron measurements for materials located in the vicinity 
of the QCEP of the FOVT are also useful for examining the possibility of whether an 
almost g-independent mode appears near q = 0. 

In this paper, we focused on Ce- and Yb-based heavy fermion systems. Even in the 
other materials, there remains a possibility that quantum valence criticality emerges, 
since quantum valence criticality arises from fluctuations of charge transfer between 
orbitals in itinerant electron systems with strong local electron correlations, which are 
quite common. Hence, it seems to have wider relevance. The experimental challenge 
to uncover new universality class in correlated electron systems is an interesting future 
issue. 
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